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Abstract Lorentz invariant violation (LIV) test is very important to study in the new physics. 
All the known astrophysical constraints either have a very small examinable parameter space, 
or are only suitable for some special theoretical models. Here we suggest that it is possible to 
detect the time-delay of ultra-high-energy cosmic-rays (UHECRs) directly. We discuss some 
difficulties in our method, including the intergalactic magnetic fields. It seems that none of 
them are crucial, hence this method could give a larger examinable parameter space and a 
stronger constraint on LIV. 

Key words: cosmic rays — gamma rays: bursts — ISM: magnetic fields — relativity 



1 INTRODUCTION 

Lorentz invariant violation (LIV) test dPavlopoulosl 1 19671) is significant for the applicability of Special 
Relativity, even if it has few trustworthy theoretical foundations. The theoretical approaches for LIV are 
mainl y from new physics, i ncluding the S tandard-Model Ext ension, noncommutative geometry dSzabol 
l2003h . loop quantum gravity dRoveliil. [l998) and string theory dGreen et all 1 1 9871: IPolcIiinskiL 1 1 998h . 

The Standard-Model Extension approach^ dColladav and Kosteleckvl l 199711 19981) is the most straight- 
forward one, which introduces LIV as an assumption. LIV may be caused by spontaneously violat- 
ing the vacuum solution, if not by the theory itself. The minimum Standard-Model Extension wishes 
to maintain all the conventional desirable properties of the Standard-Model beside allowing for viola- 
tions of Lorentz symmetry, hence it does not have many astrophysical (time-integral-kind) applications. 
How ever, the Standard-Model Ex tension can actually induce some kind of birefringence effects for pho- 
tons dKosteleckv and Mewei,l200ll) . 

Noncommutative geometry has a lot of phenomenological applications. How ever, most of them are 
based on terrestrial experiments dHinchliffe et all 120021: iKonopka a nd Maiod 20021) . The derivation of par- 



ticle Lorentz-violating terms from noncommutative geometry dCarroll et al.. 2001) seems more natural than 



other approaches, but unfortunately, it does not have (at least we don't know how it can have) a beautiful and 
feasible way to be approved by time-integral-kind experiments. The differences are as follows. Some th e- 
oretical models result in a constant space of light (e.g. by K-Minkowski space-time (Tama kTet al.ll2002l) ). 
And, researchers seem to have different opinions on whether the Lorentz-violating term 6^ depends or not 
on position, energy or momentum^. 

The loo p quantum gravity a pproach seems the most usable one. The propagational calculations 
of photons ( Alfar oet all l2002at iGambini and Pullinl 1 19991) and neutrinos (or other massive spin- 1/2 
fermions) dAlfaro et all 2000, 2002b) are fulfilled. The propagation speeds in both cases are non-trivial, 
with velocity departure linearly depending on particle energy. Further more, photons have a first order and 
neutr i nos have a second order birefringence effect. H owever, although "foam" structure (Dopl icher et all 
1 19951 iGaravi 1 19981 lHawkingi 1 19781 IWheeleil 1 1964 is really an intuitive way to understand the nature of 



* E-mail: congxin . qiu@gmail . com 

1 <http://www.physics.indiana.edurkostelec/faq.html > 

2 The kind of experiments we are interested in work only if 9^ v depends on energy and momentum, but is independent or almost 
independent of position for us to integral the effect. 
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quantum space-time, we ha ve to warn ourselves time and again that loop quantum gravity theory itself has 
some theoretical problems dAshtekar et all[l992l) . "weave" states and coarse graining approximation are at 
most effective models. 

A nother leading (an d in fact, chronologically the "first") root for the LIV calculations is from Liouville 
string dEllis et all ll992). a phenomenologic al model which makes the calculations of propag ation equation 
in the framework of string theory possible dAmelino -Camelia et al 1 1 19971 lEllis et all l2000al) . What they 
can calculate are the so-called "photons" which are the endpoints of open-strings attached to D-branes, 
and the space-time foam is described by D-brane fluctuations. The model can result in Lorentz-violating 
propagation equation by LIV of string ground state, although it also has some inconsistencies. As a result, 
LIV is stochastic, and the degree of velocity departure is first ordered. However, there is no evidence to 
support birefringence, which is in conflict with loop quantum gravity results. 

There are also some other ways to discuss LIV from the theoretical viewpoint, although some of them 
are f ormerly due to the so-called GZK anomaly, which may in fa ct be some kind of experimental er- 
rors dHiRes CollaboratiorA Eool iPierre Auger Collaboratio nl l2007l) . The methods include simply adding 
tiny (first order or second order or whatever we want) Lorentz-violating terms to a conventional Lagrangian 
(these may b e considered as some kind of Standard-Model Extension) and seeing how they can affect our 
observations dColeman and GlashowMl999tlMvers and Pospelovl 120031) . calculating th e geodesic in a topo - 
logical fluctuated classical general relativity to get some very complicated results ( Yu and Fordl [l999), 
deforming the measure of integration in Feynman grap hs (which is equivalent to inventing a new renor- 
malization skill) to g et an effective LIV dAlfaroi l2005 a.b). calculating the graviton induced corrections to 
Maxwell's equations dDalvit et alll200lT) ; however, the resultant spee d of light correction in the last method 
is independent of energy. A recent work bv lGogberashvili e t al. (2007) deduces the dispersion relation (with 
no birefringence effect) from fat brane-world scenario; but the resultant constraint seems to be too strict to 
trust that model. 

Astrophysical experimental (dis)confirmations of LIV often use far transient sources emitting high- 
energy particles. The common sources are gamma-r ay bursts (GRBs) which are cosmological, have very 
short durations, and can emit high-energy photons dGupta and Zhangl |2007|) and neutrinos 0. The other 
common sources are giant 7-ray flares of active galactic nuclei (AGNs); however, there have not been 
suitable models for the shapes of the time profiles until now. If energy can affect particle speed by the 
LIV effect (it's not the same as the effect of particle mass, which becomes unimportant if the particle is 
sufficiently energetic), as some theoretical works predicted, particles emit simultaneously from the source 
but with different energies will exhibit a time-delay w hen observed. The poss i ble ways include testing 
the time-de l ay of prompt em ission photons from GRBs dAmeli no-Camelia et al., 1998; Elli s et alll2 000b. 



20031 l2006t iNorris et all 1 1999b and giant 7-ray flare s of AGNs dBiller et al.lll999HMAGIC Collaboratioii 

2007), the time-delay of neutrinos from GRBs (lAlfaro et al. , 2000, 2002a iBertolami and Carv alho. 
ur, — r - ^ — ttt- — 1 u„„„ i rr — ; ttt — 1 , • , , . * ^7^~TT 7~n 



2000; Choubev an d Kingl 120031; [J acob and Piranl 12007b. " the polarized photons from G RBs (Fanetal. 
2007; iGleiser and Kozamehi 1200 It iG.Mitrofanovi |2003|) and distant gala x ies ( Kostelecky and Mewes, 
2001) which should be destroyed by birefringence (Alfaroetal. , 2002a; Gamb ini and Pulliri 1 19991) . 



2003) which should not be 



the synchrotron radiation from the Crab nebula dJacobson et all 12002 , 
obse rved if photons can be b oth superluminal and subluminal but electrons can only be sublumi- 
nal dMvers and Pospelovl 2003b. There are also a lot of theoretical works to explain the GZK anomaly by 
Lorentz-violating terms dAlfaro and Palmal 120031 lAloisio et al. , 2000; Ameli no-Camelia and Piranl 1200 It 



IColeman and GlashowHl999l) . so if the GZK cutoff ( Greisen , 19661 IZat sepin and Kuz'min. 1 1966 ) does in 



fact exist, the inverse proportion may also give some kind of constraints. 

The purpose of this paper is to suggest a different way to (dis)confirm the LIV effect; that is, to test 
directly the time-delay of ultra-high-energy cosmic-rays (UHECRs) from far away sources. This method 
may give a larger examinable parameter space and a stronger constraint. 



3 There are really a lot of different models for GRBs to emit ultra-high energy neutrinos, from (Waxman and Bahcall, 1997) until 
now. Nearly all of the scenarios are p + p or p + 7 7r + => v, but in different environments. See Waxman 1 2001a) for a review. 
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2 CALCULATION 

2.1 Naive Time-Delays by the LIV Effect 

One possible way to (dis)confirm LIV is simply to test the time-delay of UHECRs from far away 
sourc es. Because in mainstream q u antum gravity models, the depa r ture of velocities depends on energy lin- 
early dAlfaro et all l2000l l2002afbt |Xmelino -Camelia et al 1 119971: lEllis et all |2000a; Gambimand Pullin, 
1999) in the massless approximation, the time-delay is very sensitive to ultra-high-energy particles. A naive 
calculation shows that the time-delays are really huge. For example, in the standard cosmological model 
where Ho, Q, m and tt\ as the customary cosmological parameters, the propagation equation and time-delay 
for a massless particle has the form 



and 

A* = J_ f f^J\ (1 + *')<**' 

QG H J \ZEpJ ^n m (l + z'f + Qa 

where E <§C E v \ is the energy of the particle, z is the redshift of the source, £ is a free parameter to describe 
the degree of violation (which we want to restrict) with assumed typical value of unity, c is the speed of light, 
and E v \ is the Planck energy. To give a straightforward example, insert E — 10 198 eV ~ 6.3 x 10 19 eV as 
the GZK threshold energy, z — 0.1 ~ 400 Mpc as a nearby source, and £ = 1 as a typical dimensionless 
free parameter, we have 

At QG ~ 7yr. (3) 

We choose z = 0.1 ~ 400 Mpc rather than larger distances to avoid A£qq to be too large to be 
compared with human longevity. In this case, cosmological models are in fact irrespective, so the situation 
differs from considering less energetic but neutral particles (like photons or neutrinos) that come from more 
far away sources. Closer sources are also possible (and maybe even better), because nearly all the time- 
delay effects (including the intergalactic magnetic fields, which we discuss in detail in § 12.4b caused by 
propagation depend linearly on distance, and distance is irrespective when contrasting which one of the 
time-delay effects is more important. Remote sources are only needed when A£qq is too small compared 
with the internal duration of the events themselves, which is only several seconds for GRBs and some other 
transient sources. 

When the energy of the UHECR particles exceeds the GZK threshold, it is less possible for the source 
to be too far away, because the particles lose energy by interactin g with CMB p hotons. The main mech- 
anis ms of energy loss o n the road are photomeson production dSteckerl fl968) and e + e~ pair produc- 
tion (iBlumenthaA 11970b . with their mean free paths already being calculated. However, for the UHECR 
events with energy larger than the GZK threshold we have already observed, their time-delays by the LIV 
effect are really interesting, because they should be more energetic and more sensitive to LIV when just be 
emitted. 

However, the calculation of energy loss rate dE/dx (where x is the propagation distance) is very 
difficult, although what we have to face are trivial details of standard quantum field theory and phase s pace 



integrals. Here we simply use the existing numerical results dXh aronian and Croninl l 1 994 ICroninll 19921) to 



proceed our calculations. Time-delays depending on different propagation distances are shown in Fig. ([T] 
bottom). We see that although the time-delays finally tend to the same level as others (because all their 
energies converge to the GZK threshold energy after a long way of propagation), their differences are 
tremendous when just been emitted. So, confirming the sources of the UHECR events above the GZK 
threshold energy may also be a way to test the LIV effect, even if the source is really nearby. 

We have to emphasize here that the particle energy versus propagation distance relation in Fig. (Q] top) 
is based on some statistical results with large samples, because the photomeson interaction is stochastic 
(see more detailed discussion in § 12.4.21 ). The me an free paths fo r UHECRs with energy E > 10 21 eV are 
about 10 Mpc, and larger for less energetic ones dSteckerul968l) . A certain particle can only suffer from a 
couple of collisions before being observed, so the LIV constraint by one single UHECR event has its initial 
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Fig. 1 Top, how particle e nergies change with propagation distance, if their initial energies are 
above the GZK threshold dCronid[l992l) . The initial energies are chosen to be 10 20 , 10 21 and 
10 22 eV. Bottom, the total time-delays caused by the LIV effect. The shadow regions are the 
time-delays caused by the magnetic field of the Milky Way, which can be neglectable undoubt- 
edly in the late two cases. 



measurement errors. However, the observed time-delay can at least give an upper limit for LIV, and we can 
still improve our result by averaging more events or using more advanced statistical methods. 

2.2 Problem of the Applicability of the Time-delay Equation 

To investigate the capacity of the method we propose, there are several problems which have to be con- 
sidered carefully. The first problem is whether the naive propagation equation can be used for UHECRs. 
Although we are not sure w hat the compositions of UHECRs are, air shower data exclu de photons or 
electro ns, and the existence dHiRes Collaboration! [2008L IPierre Auger Collabor ation! 120071) of the GZK 
cutoff dGreisenl 1 19661; IZatsepin and Kuz' mini 1 1966) suggests in fact protons or heavier nuclei. Although 
we cannot rule out other possibilities like exotic particles, we will assume in the context they are protons 
which are compound and have finite rest mass (discussions with the assumption that they are heavier nu- 
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clei like Fe are analogous). Mass is not a serious problem at the energy scale of the GZK threshold^. For 
example, in the case of parameters using above in Eq. ([3), the time-delay affected by proton mass is only 
5 x 10~ 6 s, much shorter than AiQQ affect by the LIV effect. A more serious problem is the complexity of 
protons. We all know that proton is made up of three quarks, therefore using the LIV calculations for ele- 
mentary particles to calculate it will not be justified. A detailed study of quantum chromodynamics (QCD) 
with Lorentz-noninvariant terms is needed; howev er, it will certainly be very diffi cult. We still use the same 
propagation equation by some scaling arguments dColeman and Glashowi 1 19991) . or by simply thinking its 
effect as an overall constant coefficient, just like that of turning off QCD0. For the reason that the LIV 
confirmation is still qualitatively rather than quantitatively at present, an overall constant can be neglected. 



2.3 Problem of the Source 



2.3.1 Time Bases 



The second problem is how we can choose the time bases, thus compute the time-delays with a suitable 
zero point. Because At should be typically very long, as shown in Eq. (01, it is a serious problem to know 
when UHECRs should come if we change their energies (because A£qq depends on energy of the particle) 
or turn off the LIV effect (as in the classical limit of E — > 0). Comparisons of events with photons or 
other low energy massless particles (which are certainly much less energetic and can be taken for as in the 
E — ► limit; hence the LIV effect of it can be neglected) and other UHECRs (with energy different from 
each other) from the same source can scale the time-delay. The precondition is that we have confirmed the 
source, or that we are convinced of t he fact that different signals come f rom the same source. 

As the recent powerful evidence dPierre Auger Co llaboration!, 120071) shows, UHECRs come from some 
extragalactic sources because they are anisotropic and correlated with the direction of the Super-galactic 



plane. In that case, the mainstream models for sources are AGNs (Ginzburg and Svrovatskii, ll964tlHillasL 
11984 and GRBs dVietril 1 19951; IVietri et all 120031: IWaxmanl 1 19951 l2004t IWick et all |2004 . but other 
sources distributing within the Super-galactic plane are also possible (if they are related to, e.g. galaxy 
formation or stellar formation, which is always true). The main mechanism is Fermi acceleration but in 
different environments. 

As AGNs are lasting sources, we can hardly know very well when the source emitted the UHECR 
particle we observed. However, some recent theoretical work dFarrar and Gruzinovl 120081) shows that 
the U HECR emissions ar e associated with AGN giant flares, with typical wait-time about 10 3 to 10 4 
years dDonlev et all 120021) . Because the duration is much longer than the typical time-delay we gave in 
Eq. (0), AGNs can be used in this method if the theoretical work mentioned above is true. GRBs are much 
better sources, because nearly all of the mainstream central engine models (including collapsars, supranovae 
and mergers of compact objects) tell us they are transient and burst only one time in their whole lives. If 
the emission of UHECRs and the burst itself happen almost at the same time (which is the most reasonable 
assumption), we can scale the time-delay by the observed low energy 7-rays because they can hardly be 
affected by mass, electromagnetic fields and the LIV effect. Other sources are also possible, if they emit 
particles (photons for instance; however, they are not exclusive choices) other than UHECRs which can be 
observed by our scientific equipments. 

Are these kind of sources practicable for our purpose? The distances of the sources mentioned above 
are all suitable for the constraint that A£qg given in § 12.11 should not be too large. Short GRBs are often 
not too far away from us, and there are already a number of GRBs wit h redshift z ~ 0.1, including a 
special one (GRB 980425) with an especially small redshift z = 0.0085 dGalama et all fl998l) . Although 
the number density of AGNs decreases quickly when z < 1, there are already hund reds of nearby AGNs 
have been observed until now (e.g. the V-C catalog dVeron-Cettv and Veronl [2006b has 694 AGNs with 
redshift z < 0.024). Similarly, it is reasonable to assume that other possible sources of UHECRs are not 
too far away, because UHECRs are not isotropically distributed in the celestial sphere. 



4 Of course, massive and massless particles are totally different from the viewpoint of quantum field theory. However, we give up 
discussing it deeply, because of the still inconsistent theoretical works. 

5 If we turn off QCD, the propagation equation can be used to every elementary particle inside a proton, so the overall effect is 
only a constant coefficient. 
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2.3.2 Confirmation of the UHECR Sources 

The assumption in the above paragraph is that we have confirmed the source, or we are convinced of the 
fact that different signals come from the same source. However, it is not always the case. Notice the fact 
that UHECRs are singular events (it seldom happens that the UHECR events have clustering properties), 
confirm their sources by statistical correlation is very important. 

Metrical bias in spatial dimensions are caused by (i) intergalactic magnetic fields and (ii) the uncertain- 
ties of detectors; the LIV effect cannot affect the orientation of UHECRs. If the collective effect of (i) and 
(ii) is small enough, we can confirm the sources by their locations in the celestial sphere; however, it may 
not be the case. If we assume that the effects of (i) and (ii) are both stochastic, confirmation of the sources 
is a pure statistical inferential problem. Astrophysical parameters only affect the statistical sample s by (i) 
the UHECR energy band or (ii) the possible correlative time interval. iPierre Auger Collaboration! (120071) 
has already discussed the statistical correlation between the arrival directions and the positions of known 
AGN. The same method can be used for our purpose; however, their arguments do not include the temporal 
dimension. When discussing the LIV effect, temporal dimension is very important. Hence, we should put 
by hand a possible correlative time interval when choosing the statistical samples; that is, assume that the 
collective time-delay caused by LIV, intergalactic magnetic fields and other reasons does not exceed this 
interval. Notice the fact that the observational history of UHECRs and correlative sources are at most sev- 
eral decades, which may be shorter than the collective time-delay, it is a good idea to ignore the temporal 
dimension and choose all the samples we know to do the statistical correlation. However, if the intergalactic 
magnetic fields are sufficiently large, we will never know the sources of UHECRs, no matter whether LIV 
exists or not. 



2.4 Problem of the Intergalactic Magnetic Fields 

The third but the most annoying problem is the intergalactic magnetic fields. Because protons take charges, 
their trajectories will be (Larmour) curved by magnetic fields, and the departures from straight lines will 
cause extra time-delays. Our method is only suitable when the time-delay Atu by the magnetic fields is 
less than by the LIV effect. 

Because an UHECR particle keeps constant energy inside some homogeneous magnetic field, the time- 
delay should be 

1 D 3 ( D \ 3 / E \ ~ 2 / B \ 2 

where D is the linear distance of the trajectory, B± is the perpendicular magnitude of the magnetic field, E 
is the energy of the particle, and = E/(c ■ eB±) is the Larmour radius. 



2.4.1 Comparison with Photons 

For simplicity, we first discuss the way of comparing the UHECRs' time-delay with photons, because 
photons are irrespective of the magnetic field, and their time-delay by the LIV effect can be neglected 
compared to UHECRs for their relatively lower energies. 

The real trajectory can be devided into three parts, inside the host galaxy, inside our Galaxy and in the 
intergalactic media (IGM), that is, A^m = AiM.host + A£m, Milky + AiM.iGM- We have already chosen the 
values of D and B± both for a typical galaxy in Eq. (0), so AiM, Milky ~ 0.79 yr is the typical value for the 
time-delay effect of the Milky Way, which can be negligible compared to AiQQ we estimated in Eq. ©. 
Of course, A£qg decreases when the source comes nearer, but the effect by the Milky Way's magnetic 
field remains unaltered, so it would be troublesome when considering the use of more nearby sources to 
test LIV, as mentioned in S I2.1I However, because the time-delay by magnetic fields is absolutely classical, 
when we have fine structure models for the magnetic field of our Galaxy someday, we can deduct this effect 
directly @. When the UHECR particles are initially more energetic than the GZK threshold, effect from the 



6 Because the correlation length of the magnetic field in our galaxy should be compared with the scale of the galaxy itself. Further 
more, we know the direction where the UHECR particle is related to the local structure. 
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Milky Way's magnetic field can always be negligible, as shown in Fig. ([T] bottom). The time-delay by the 
host galaxy of the GRB will not be worse than by our Galaxy, because the energy E will be larger (if it 
formally exceeds the GZK threshold) or at less equal (if less than the GZK threshold) when just emitted. 

However, the effect by the large from scale intergalactic magnetic fields is more thorny, because until 
now we lack goo d models for the ma gnitude and topological structure of the fields. A constraint from the 
CMB anisotropy dBarrow et all 1 1997b gives 

Bi GM < 6.8 x I0- 9 (n h 2 ) 1/2 G - 4.9 x 1CT 9 G, (5) 

where we choose fio = 1 and Hubble co nstant Hp = 72 kms~ 1 Mpc -1 . Another constraint from the 
observed rotation measure (RM) of quasars (Rronbergl ll994l) gives 

B lGM < 10- ( T A^) G, (6) 

where A denotes the correlation length (coherence length) of the magnetic fields, as a reasonable assumption 
that the power spectrum of magnetic fields has a large scale cut-off. 

If we assume that the field is conglomerated and homogeneous inside every segment (with typical 
scale of correlation length A), the UHECR particle will randomly change its direction due to Larmour's 
motion, but goes a nearly strict line as a whole. If £>igm is independent the correlation length A, the overall 
time-delay should be 

D \ ( A \ 2 / E \ 2 f B\ ' 2 



When B depends on the correlation length as B — BqX^ 1 ^ 2 G in Eq. ©, the time-delay is 

Ai M ,iGM * 1.18 (40^) (7^) ( 6 . 3x f i9 eV ) (iFTTc) < 8 > 

Other possible parameters are denoted in Fig. (0. We see that it is neededQ for our purpose that £?igm is 
slightly less than the upper limits given by Eq. (|5]l and (O, unless we choose a smaller correlation length. 

Things will be more interesting when consider the UHECR particles with energy exceeding the GZK 
threshold. Noting that B± (£>o.±) and A in Eq. (0 and (O are independent of the source properties, we may 
define 

— Vf-^ ' 

IMpc/ vio _11 g 



^1-18(tT7Tt) [77^77) W 



^ 1-18 X ( VttSt77 ) <10> 



in Eq. (0 and 

n = 1 18 y I 

IMpcy ViO -11 G 

in Eq. ([8]); then the effect by intergalactic magnetic field has a uniform expression 

A£m,igm ^ V ■ (45^) ( 6 . 3x f i9 eV ) y r - (») 

In Fig. (0, we calculated A£qg + A^m, Milky + AtM.iGM in all, with r\ = 1, 70 and 5000 respectively. 
T) = 5000 has already saturated the upper bound given by Eq. (0 and ©, so A^m.igm cannot be larger. 
Noting that when E > 10 21 eV, the UHECR particle will absolutely not be affected by magnetic fields if 
it is not too far away (roughly D < 10 Mpc), hence the only thing that can make a visible time-delay is 
the LIV effe ct. HiRes and AGASA have already observed a couple of the UHECR events with energy E > 
3 x 10 20 eV (lAGASA Collaborationil200llHiRes Collaboration[|2005l) . If their distance D > 20 Mpc, their 
initial energy Eq will exceed 10 21 eV, as shown in Fig. ([T] top). Hence, seeking the sources of UHECRs 
with energy E > 3 x 10 20 eV will tremendously help us to (dis)confirm the LIV effect. 



It is possible in principle because Bjqm remains largely unknown by the intrinsic observational difficulties i Beck et al., 1996). 
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Fig. 2 Time-delay AiM.iGM caused by the intergalactic magnetic fields when z = 0.1 ~ 
400Mpc and E = 6.3 x 10 19 eV, with different correlation length and strength of the mag- 
netic fields. The horizontal line is A£qq = 7yr, as the example shows in Eq. (01. The method is 
useful when Ai M ,iGM < Ai M < A^qq. 



Things will be worse if the re exists a global cosmic magnetic field, or the fields have structures like 
filaments or sheets (lRv^.ll998h . Those will cause larger time-delays. The first trouble can be easily seen 
from the expression of A£m,igm given above, because it is equivalent to a huge A. The second trouble is 
because, when a magnetic cloud collapses to 1/k of its diameter, the magnetic field strength will increase 
k 2 times than it used to be. Although the particle will miss a lot of clouds it used to bang into (only for the 
case of filaments but not sheets, which make things worse), we still have 

B l x2 -» v(k 2 B ± ) 2 ( ~\) = kB 2 \ 2 , (12) 
k \k J 

and the same for the B\ A cases. 

However, in fact the irrefutable observed anisotropy of UHECRs (Pierre Auger Collaborat ion! 120071) 
has already given us an upper limit for B±_ and A, and also whether the magnetic field has already collapsed 
to filaments or sheets or not. Note that 20 among 28 highest energy events detected by the Pierre Auger 
Observatory are within a 3.1° circle of nearby AGNs (with distance less than 75Mpc away). No matter 
whether we believe these UHECR particles to originate from those AGNs or not, it is unassailable that 
UHECRs are unisotropically distributed, and seem to be correlated to the Super-galactic plane. So, the 
angular dispersion caused by the intergalactic magnetic fields should be less than a couple of degrees. The 
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Propagating Distance (Mpc) 

Fig. 3 The thick solid lines and shadow regions are the same as in Fig. (HJ. In addition, we 
calculated Aiqc + A^m, Milky + ^MiGM with rj = 1,70 and 5000 respectively. It seems that 
finding the sources of the UHECR events with energy J?> 3 x 10 20 eV will tremendously help 
us to (dis)confirm the LIV effect. See the context for detail. 



angle departure inside some homogeneous field is 

a~-^-, (13) 

and the different irrelevant magnetic bulks (with typical size A) can be considered as a random walking 
process. The overall angle departure is 



A ID _ \/D^\ 

Choosing D = 100 Mpc as the typical scale of Super-galactic plane, we have 

D \ ( E \ _1 / B ± 



(14) 



^j^ 3xiol9eV j ( vIF ^j, (15) 
and 

aiGM " 4 ' 21 °(ioo^) im^) (e.sxfo^ev) (iF^g)" (16) 

Notice that when oitgm approaches a couple of degrees, as the above equation shows, the upper limits given 
by Eq. (O and © have already been saturated. In addition, because of the fact that A 1 / 2 ^^ — > VkX 1 ^ 2 B±, 
filaments or sheets can also be suppressed. 

One question is whether the Pierre Auger data tell us that oigm should be equal to (rather than less 
than) several degrees? Absolutely not. oiigm can a ls° be much less (so B and A can also be much less). 
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Even if we have confirmed some UHECR sources, the angular dispersion can also be caused by reasons 
other than o;igm> for example, magnetic field of the Milky way or simply the measurement errors. 

We notice that some authors g ave a larger A£m,igm compared to ours given in Eq. © and ©. 
IWaxman and Miralda-Escude1(ll996l ) gave AtM.iGM ~ 100 yr becaus e their corre lation length A ~ lOMpc 
is 10 times larger than ours (equivalent to T) = 100 in our definition). ISigll (l200ll) gave A<m,igm ~ 10 3 yr 
because he chose a re ally large magneti c field B ~ 10 9 G (equivalent to r\ = 10 4 ); however, with a smaller 
travelling distance D. IWaxmanl (1200 lbl) gave an upper bound of AtjvijGM even as large as 10 7 yr, because 
his typical magnetic field B ~ 10 s G is really huge. He also argued that A£m,igm > 100 Y r by some 
statistical reasons of nearby source candidates and the UHECR events above the GZK threshold. The first 
two estimations are consistent with our constraint from correlation of Super-galactic plane and the UHECR 
events, the few discrepancies are just because we choose different typical parameters (which are all possi- 
ble according to our current knowledge, because we know really little about the true value of B and A) to 
write our formulas. We suggest that the anisotropy of UHECRs can give a tighter constraint of intergalactic 
magnetic field strength B, so the upper bound of A^igm we can assure at present should be as low as 
10 4-5 yr. The lower bound A^m,igm > 100 yr can be overcome because we know really little about both 
possible nearby sources and the UHECRs events, and the estimation is dependent on some details of source 
models. In addition, all the estimations given above are only suitable for particles with energy below the 
GZK threshold, because the energy loss is ignored. As we show in Fig. (0, the effect of intergalactic mag- 
netic field is much less important if the energy of the observed UHECR event is much larger than the GZK 
threshold. 

2.4.2 Comparison with Other UHECR Events 

We can also compare the time-delay with other UHECR events (with slightly different energies), emitted 
nearly simultaneously from the same source. Of course, because the UHECR events are really rare, it may 
hardly happen. 

In this case, blurs in both arrival direction and time-delay have to be analyzed carefully, (i) Blurs have 
two reasons. Particles with different energy follow different trajectories, thus leading to different directions 
and time-delays, because of the random topological distributions of the intergalactic magnetic fields, (ii) 
At the same time, particles above the GZK thre shold energy would interact with CM B photons as the 
Poisson processes, introducing extra randomicitv. I Waxman and Miralda-Escu de (1996) discussed the blur 
effect with UHECRs below the GZK threshold, in which case energy loss by photomeson production can 
be ignored. At the end of § 12.11 we have already discussed a little the influence of the LIV time-delay by 
stochastic photomeson production. 

It is easy to understand that when the particles are extremely energetic, blurs in both arrival direction 
and time-delay caused by the intergalactic magnetic fields become less important. However, using one 
of the UHECRs to scale the others may be dangerous, if their energies are large enough (e.g. larger or 
equal to the GZK threshold) to make us believe that they have suffered one or more times the photomeson 
interactions. Because of the randomicity of Poisson arrival photomeson interactions, particles observed 
with the same energy from the same source may have tremendously different interacting histories and thus 
possess different time-delays caused by both the intergalactic magnetic fields and the LIV effects. 

However, for UHECRs less energetic than the GZK threshold, photomeson production is turned off, 
and comparison becomes possible. The requirement that the intergalactic magnetic fields should not be very 
large, is the same as in the case of comparing UHECRs with photons, which we have already discussed in 

§E33] 

2.5 Problem of the Energy Measurements in Air Shower Detectors 

Notice that the energy measurements in different mass composition of the UHECR events and different 
air shower detectors have disagreements from each other which cannot be negligible, so it is necessary to 
discuss here the influence of the LIV confirmation by energy demarcation uncertainties. In § 12.41 we have 
already discussed two different methods to restrict LIV, the comparison (i) with photons and (ii) with the 
different UHECRs respectively from the same source. 
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For the reason that the investigations of the LIV confirmation are qualitatively rather than quantitatively 
at present, the energy metrical uncertainties are not crucial for method (i), because it can only introduce an 
order one coefficient of £ in Eq. (fJJ. When the UHECR events are not too energetic to neglect the time- 
delay caused by the intergalactic magnetic fields, absolute energy measurements are important. However, a 
global constraint for the collective influence of Aiqc + A£m,igMi hence the upper limits for both A£qq 
and AiM,iGM respectively, are still suitable for our purpose. 

For method (ii), things are a little more complicated. Uncertainties introduced by the different assump- 
tions of mass composition are not crucial. The reason is that, when assuming different UHECRs to be the 
same kind of particles (protons in our context), a mistaken assumptive mass composition can only introduce 
an order one coefficient of £, just as in the case of method (i). However, it is intractable for UHECRs de- 
tected by different air shower detectors with different energy metrical techniques. A wiser way is to choose 
some kind of calorimetri c measurements to determine UHE CRs' energies by different detectors (like fluo- 
rescence light emissions dLinslevl Il983t ISong et all 120001) ) which are relatively model independent. As a 
matter of fact the discussions of LIV are presently still superficial, we may hope that energy demarcations 
are finer for further investigations of LIV in the near future. 

3 DISCUSSION 

3.1 Two Known Events 

There was an archaeological report about the association of UHECRs and GRBs (Milgrom a nd Usovlll995l) . 
The authors found that GRB 910503 and 921230 are associated with two highest-energy cosmic -ray shower 
events, with really small error boxes and time-delays of 5.5 and 1 1 months respectively. If GRBs are really 
sources for those two UHECR events, there are very strong constraints both for LIV and the strength of in- 
tergalactic magnetic fields (as the time-delay is much shorter than the naive estimation we make in Eq. (01), 
because all effects such as rest mass, magnetic fields and quantum gravity, are addible, and to ignore some 
of them gives the upper constraint for the rest ones. However, we should not be too serious for that kind of 
stories, because they may only be a coincidence. 

3.2 Comparison with Other Models 

Although there are other constraints of LIV which are much stron ger than the method we suggested, the 
method mentioned above h as also its special purpose. Birefringence dFan et alll2007t lGl eiser and Kozameh. 
2001; G .Mitro fanovl 120031) can only be calculated in the framework of loop quantum gr avity but not in 
other appro aches, hence it may be wrong in a whole. The synchrotr on radiative constraint dJacobson et all 
I2002ll2003l) depends on a special theory dMvers and Pospelovll2003l) . which needs a dimension-5 Lorentz- 
violating terms to induce birefringent photons but subluminal electrons (whose maximum speed cannot 
converge to c). The inverse proportion of the GZK anomaly may also give some stronger constraint. 
However, the scattering dynamical discussions are always only one-sided, which means that a scattering 
channel is open or suppressed only if the effect of LIV in opposite for two relative particles (and therefore 
their velocity as well as their effective mass being different). Although in the old days, the GZK anomaly 
is the most important reason for theoreticians to study LIV, its inexistence (IHiRes Collaboration! [2008; 
iPierre Auger Co llaboration. 2007) has not borne down the LIV subjects. 

Testing the time-delay of UHECRs is a more direct way to study LIV. It can contain most kinds of the- 
oretical works. If the intergalactic magnetic fields are sufficiently small (which is still absolutely consistent 
with the observations until now), it may have larger examinable parameter space for violation scale £ (in 
Eq. (O) than using photons or neutrinos. Even if its examinable parameter space is in fact much smaller, 
for the reasons mentioned above, the other causations are all classical and thus can someday be subtracted 
by models. 

4 CONCLUSIONS 

We have suggested to (dis)confirm LIV by simply detecting the time-delay of UHECRs. We considered 
some other reasons which also cause the time-delay, including the intergalactic magnetic fields. If the en- 
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ergy of the UHECR events we observed is below the GZK threshold 6.3 x 10 9 eV, a typical intergalactic 
magnetic field B < 1CT 11 G and correlation length A < 1 Mpc may be needed to give a parameter space 
examinable enough to constrain LIV. However, for an UHECR event with energy larger than 3 x 10 20 eV, 
our method is always possible. Because of the fact that we know really little about the intergalactic mag- 
netic field's strength, if it is much smaller than the current upper limit B < 1CT 9 G, our method may give 
a larger examinable parameter space and a stronger constraint of LIV than other constraints. 
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